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SUMMAKY 

ran a ^??? er K 11z ^ d method o 1o 9y to structural life prediction, design, and 

f a<? ! d u 5°? 3 fa ? l9ue cr ^er1on Is advanced. The life prediction 
methodology Is based In part on work of W. Welbull and G. Lundberg and A 
Palmgren. The approach Incorporates the computed life of elemental stress 
volumes of a complex machine element to predict system life. The results of 
coupon atlgue testing can be Incorporated Into the analysis allowing fo? ?1fe 

t1ca] C certa1nty C0mp0nent ° r structural renewal rates with reasonable statls- 


INTRODUCTION 

Design of machine elements Is based for the most part on yield stresses 
ertils 11 ^ 6 11m1t1n9 stresses. In addition to the components material prop- 
erties, proper consideration must be given to the effects of notches surface 

ronmlnta i * .^° l ! lponen t size, residual stress, temperature, duty cycle and envl- 
onmental factors such as corrosive or chemical exposure. For most machine 
e.ements, Individuals and organizations usually develop design methodology 
based upon engineering fundamentals found In most machine design text* and 
factors based upcn their corporate experience and test data. As a result It 
Is not too unusual fir different organizations or Individuals, starting with 
sai " e Mi Slm *4 d ? s1gn redu1reme " ts to reach dissimilar conclusions or 
to S the S prob1em Seem1n9 y app1y1ng the same fundamental engineering principals 

Setting aside the subjective creative aspects of design, there aDDear? tn 
be nonuniformity of data from which numbers and design factors are selected as 
ZllV d1fferencas the computer codes and boundary conditions used 5n the 

to i?tahmh S r; + /° C ?? p ?i! nd these difficulties, fatigue data which are used 
to establish fetlgue limits are usually of a limited nature with the conditions 
under which the data were obtained not adequately defined or reported. Such 
temperature, humidity, number of specimens, specimen size and volume 
heat treatment, hardness, surface finish and life distribution are not glOen 

The established fatigue limit for much of the reported data Is a mean 
value. From a statistical viewpoint, the median value Is equal to or less than 
rfllf ?!f'j IJ 1s c ? n be 1r, terpreted as meaning that before a fatigue limit Is 
th fh 6 a P roba b1Hty that 50 percent or more of the specimens had 
failed. In other words, even at the fatigue limiting stress, life Is finite. 

Of course, experienced design engineers have recognized this for years Thev 
factors „ t0 the1r design procedure usually based oS eipeJ ^ 

Hp^nnoH Wh l 6 t f es 5 P roced ures are generally adequate they can result In over- 
designed, oversized, overweight and overcost structures. 
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A fundamental principle of good design Is to recognize that any structure 
can fall. However, should It fall. It should be designed not to cause personal 
Injury or secondary damage. In other words, If the structure falls, It falls 
In a safe or benign manner. Once a structure Is designed to fall In a benVjn 
manner, It then can be designed for finite life whereby the overall size, 
weight, and cost can be reduced and still meet the reliability requirements of 
the application. 

In view of the aforementlon. It becomes the objective of the work reported 
herein to advance a generalized methodology for structural life prediction, 
design, and reliability based upon a fatigue criterion. The life prediction 
methodology Is based In part on the work of W. Welbull (refs. 1 and 2) and G. 
Lundberg and A. Palmgren (refs. 3 to b) . 


SYMBOLS 


C dynamic load capacity 

c stress-life exponent 

e Welbull slope, 1/a 

F probability of failure 

f failure probability density function 
h exponent 

L life 

L a adjusted life at a 90-percent probability of 

survival 

l. m mission life 

L-| o 10-percent life, life at which 90 percent of a population will survive 

n load life exponent 

P applied load 

P(x) probability for occurrence of extent x 
R cumulative renewal function 

r renewal density function 

S probability of survival 

S] system probability of survival 

1 dummy variable of Integration 

t time or time function 
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V stressed volume 

X time function or stress function 

X|j characteristic strength 

X v stress below which no specimens fall 
z depth to critical maximum shearing stress 

<x shape parameter 

0 characteristic life 

Y mean time to failure 

n stress cycles to failure 

u location parameter or time below which no failure Is expected to occur 
o u fatigue limit 

t critical stress 


STATISTICAL METHOD 

In the late 1930's (circa 1937) W. Welbull In Sweden attempted to graph- 
ically linearize various types of statistical data distributions for small 
sample sizes. By trial and error, Welbull found that by having the lnln 1/S 
as the ordinate and In X as the abscissa, where S Is the probability of 
survival or statistical percent of samples surviving and X Is a time func- 
tion, most engineering data distributions will plot on a straight line 
(ref. 7). Hence, It became possible fnr small amounts of data to estimate a 
generalized population distribution for a population of Infinite size. Having 
empirically determined this, Welbull developed a theoretical basis for what 
was to become known as the "Welbull distribution" or "Welbull plot" which was 
published In 1939 (ref. 1). Welbull (ref. 8) defines the distribution function 
as "an adequate expression for a large class of phenomena which have the prop- 
erty that the probability of nonoccurrence of an event Is equal to the product 
of the elementary probabilities," 

S . exp|- [(X - u)/0] 1/a | (1) 

The function Involves three parameters, 

a the shape parameter 

0 the scale parameter 

v the location parameter 

Where a the shape parameter Is equal to 1. 0.5 and 0.28, the respective 
distributions approximated are exponential, Rayleigh, and normal (Gaussian). 

A typical Welbull plot for rolllngelement bearing fatigue data Is shown In 
figure 1. The slope of the line "e" which Is called the "Welbull slope" Is 
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equal to 1/a. For most rolling-element bearing data e equals 1.1. The 
location parameter m Is a finite time under which there would be a zero 
probability of failure to occur. If failure can occur In one stress cycle 
then the location parameter v would be assumed to be zero. 

For fatigue data, If n stress cycles to failure Is substituted for X 
in equation (1) and e for 1/a, and where F * 1 - S, equation (1) can be 



The characteristic life 6 Is the 63.2 percent failure life of the pop- 
ulation distribution. 

Welbull stated that the probability of survival S could be expressed as 

in 1 c e u 

in s ~ t n V ( 3 ) 

where V Is the volume representation of stress concentration referred to 
herein as stressed volume, t Is the critical shearing stress and c Is an 
exponent denoting a stress-life relation where 

n ~ t- c ( 4 ) 

The values of c and e can be determined experimentally. 

The effect of stressed volume can be Illustrated wherein two specimens of 
stress volumes V-| and V 2 respectively are subjected to equal stress 
t. If ni Is determined at a probability of survival Si, then the 
probability of survival S 2 for V 2 for the life Is given by the fol- 
lowing expression 

V /V 
2 1 

S = S ( 5 ) 

2 1 


where V 2 > Vq . 

For the same probability of survival, the specimens with the larger 
stressed volume will have the lower life. This principle was applied to suc- 
cessfully normalize rolling-element fatigue data by Carter (ref. 9) and 
Zaretsky (ref. 10). This principle was applied by Grlsaffe (ref. 11) in a 
Welbull analysis of shear bond strength of plasma-sprayed alumina coatings on 
stainless steel (fig. 2). Grlsaffe showed that the calculated mean bond 
strength decreased with Increasing test area In accordance with equation (5) 
(fig. 3). Further, the stress at which no specimens failed could be determined 
by restating equation (1) as follows: 
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( 6 ) 


F x ■ 1 - exp 



e 


where the 
Fx 


X 


X v 


parameters are now defined as 

statistical percent of specimens which when tested at one sot of 
conditions, fall at given stress or lower 

stress 

stress below which no specimens failed 
characteristic strength 
Wei bull slope 


LUNDBERG - PALMGREN THEORY 

Lundberg and Palmgren were contemporaries of Welbull. A problem existed 
In the rolling-element bearing Industry to establish the lives of these bear- 
ings short of extensive testing. Taking the Welbull approach a step further, 
Lundberg and Palmgren (ref. 3) took equation (3) and added another element, 
the depth to the critical maximum shearing stress, z where 



c e u 
t n V 

„h “ 


( 7 ) 


The rationale for Including the depth z was that the Initiation of a 
fatigue crack occurred at the depth z and that the distance the crack needs 
to travel to the surface until a spall (pit) occurs Is equivalent to a critical 
crack length. In other words the greater the distance the maximum shearing 
stress below the surface the longer It takes for a fatigue crack to propagate 
to the surface, the longer the fatigue life. Lundberg and Palmgren took the 
rolling bearing geometry, kinematic, stress theory of Hertz (ref. 12) and 
Thomas and Hoersch (ref. 13) and Incorporated then Into equation (7) to obtain 
a series of equations which relate the rolling-element fatigue life of a bear- 
ing to the applied load wherein the resultant life would be In millions of 
Inner race revolutions. Knowing the speed of the bearing the life In hours 
can be determined. Lundberg and Palmgren further refined their approach to 
bearing life prediction to Include a fictitious load designated C, the 
dynamic load capacity. The dynamic load Is the theoretical load which when 
applied to the bearing would result In a life of one million Inner-race revo- 
lutions and where 


C = P 



( 8 ) 
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P applied load 

Lio 10-percent life, life at a 90-percent probability of survival, 
Inner-race revolutions 

n load-life exponent usually taken as 3. 

Equation (8) can be written: 

C n 

ho - <P> 


( 9 ) 


Hence, the predicted life of a bearing can be determined by knowing C and P. 

The Lundberg-Palmgren theory Is now an International standard used by 
every manufacturer of bearings around the world (ref. 14). The theory has 
been applied to other machine elements and mechanical transmission systems 
(refs. 15 to 28). Table 1 Is a list of applications of the Lundberg-Palmuren 
theory to life prediction. 


SYSTEM LIFE AND RELIABILITY 

The life and reliability of a system Is based upon the lives and reli- 
abilities of all of Its components. Having determined the life of the Indi- 
vidual components using Lunberg-Palmgren, the probability of survival of the 
entire system Is as follows 


S T * S 1 • S 2 • S 3 ... S n (10) 

The system life-reliability equation can be written as follows 



where L-j , L 2 , L 3 . . . L n are the lives of each component of the 
system at a 90 percent probability of survival. The system life L can be 
determined at each system probability of survival Sj. For a 90-percent 
probability of survival 


. ( 1 


\ 2 (\ 

e o e 

\ 3 /| \ n 


1 *0 

;) *0 


;) -•(y 

(12) 


MISSION LIFE 

A system does not usually operate at one constant load In actual service. 
Miner's Rule Is used to sum fatigue damage of a mission profile consisting of 
loads and tlme-at-loads. For a given probability of survival, the mission 
life for the system L m Is 


L 


m 







-1 
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(13) 



where t a , th, t c and tj are the fraction of the total time at loads 
Pa. Pb. Pc a nd Pd. respectively, and L a » Lb, L c and Ld are the system 
lives at a given probability of survival at loads P a , Pb. Pc and Pd, 
respectively. Welbull plot or distribution can be constructed using this 
method by determining the mission life at varying probabilities of survival. 
This method was used by Lewlckl to determine the fatigue life of a turboprop 
reduction gearbox (ref. 28). The analysis was compared to actual field data. 
The results of this comparison are shown In figures 4 and 5. 


COMPONENT REPLACEMENT RATE 

The Welbull analysis Is a valuable tool for predicting the life of com- 
ponents or systems. The analysis would describe the failure rate In field 
service only If all components were put Into service at the same time and If 
failed components were not replaced. But, a certain number of components must 
be kept In operation, failure rates with replacements Is of Interest. As an 
example, If there are 10 000 bearings In the field, the user or manufacturer 
must know how many spare bearings will be needed over a given period In order 
to keep 10 000 bearings In operation at all times. Over the Interval of ser- 
vice operation It may take 15 000 replacements to keep all 10 000 bearings 
running. The Welbull analysis can never exceed 100 percent, but field failures 
can and often do exceed 100 percent (ref. 19). 

The cumulative probability of failure for the first time, assuming con- 
stant service conditions, Is a function of the length of time that the bearing 
has been running. Failure Is the complimentary function to survival according 
to the following relation: 


F(t) = 1 - S(t) (14) 

The probability density function for failure Is 

«*>■£ < 15 > 

The Instantaneous probability of failure for the first time In any time Inter- 
val from t to t + At Is then given by 

P(f1rst failure In At Interval) = f(t)At (16) 

By using renewal theory (refs. 29 and 30), the probability of having to replace 
a bearing In the field Is written as follows: 

P(mak1ng a replacement) = r(t)At (17) 

where r(t) Is called the renewal density. The renewal density Is calculated 
from the following summation (ref. 29): 

X 

r(t) = £ f (t) (18) 

k=l K 

where f|<(t) Is the k-fold convolution of f with Itself and Is computed 
by the following recurrence relation (ref. 29): 
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W t3 


(f, 5 f) 


(19) 



- T)dT 


The expression f k( t) At gives the probability of a k th failure occurring 
In the time Interval t to t + At. Since, when a failure occurs, 

It can be for any value of k, It follows that the renewal density function 
should be defined as the sum over all the f^'s. 

The total number of replacements made during the first t units of time 
Is obtain by Integrating the renewal density function as follows: 

R(t) r(T)dT (20) 


For a group of bearings that has been operating for some time with failures 
occurring and replacements being made, It Is also Important to know the MTBF 
(mean time expected between failures). According to the renewal density 
theorem (ref. 30), 

11 m r(t) - 3 = ; (21) 

t-« jf* tf (t)dt y 

Therefore, the mean time between failures Is calculated as 


MTBF 


total number of bearings 


( 22 ) 


Reference 19 reports the use of a computer program to evaluate 
equation (18) for rolling-element bearings. A Wei bull slope of 1.1 was 
assumed. Figure 6 shows plotted results, which give the renewal density func- 
tion for a case of where the failed bearings are removed from service and 
replaced with new or restored bearings. For comparison, the probability den- 
s', .y for failures with no replacement (Welbull density function) Is plotted 
also. The area under the curves represents the probability of failure. 


The functions plotted In figure 6 were numerically Integrated and are 
shown plotted In figure 7. These are the cumulative functions for renewal or 
failure. The cumulative renewal functions Indicate 100 percent replacement 
bearings will have been needed by the time 8 L-|g Intervals have elapsed. 

By comparison, at the time of 8 Liq the Welbull distribution shows only 
65 percent of the original bearings will have failed. The difference In total 
failures would be due to replacement bearings falling. 

The mean time between failure (MTBF) from equation 22 can be obtained 
from figure 6 . The MTBF Is the Inverse of the probability density for failure. 
As an example a probability density for failure of 0.12 would give a MTBF value 
of 8.3 L-) 0 . 


FAGIGUE LIFE MODELING 

The Welbull and Lundberg-Palmgren analyses have been primarily applied to 
high cycle fatigue with the material subjected to a Hertzian stress field. 
However, as Indicated In reference 11, the Welbull analysis can be applied to 
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other types of durability problems. What Is Important Is that a material 
element only knows the state of stress It Is subjected to and not whether It 
Is In a bearing, gear, shaft, compressor or turbine. The crack propagation 
time In high cycle fatigue Is generally a small fraction of the total time to 
failure. In low cycle fatigue, the crack propagation time Is generally a sig- 
nificant fraction of the time to failure. In either case, the end result Is 
total fracture of the component making It no longer useful for Its Intended 
purpose. If the failure distribution Is within standard statistical ranges, 
then It can be represented by the Welbull analysis. Hence, the analysis should 
be blind to high or low cycle fatigue. This was recognized by loannldes and 
Harris (ref. 31). Using Welbull and Lundberg- Pa Imgren, they Introduced a 
"fatigue limiting stress" and Integrated the computed life of elemental stress 
volumes to predict life. This approach leads to a method of applying Lundberg- 
Palmgren life prediction techniques to other components besides those subject 
to Hertzian loading. It also allows an Investigator to use the results of 
coupon testing to predict the life of complex shaped components subjected to 
nonHertzlan cyclic stressing. loannldes and Harris (ref. 31) applied their 
analysis successfully to rotating beam fatigue (fig. 8), flat beams In reversed 
bending (fig. 9) and beams In reverse torsion. Based upon their approach or a 
modification thereof, design procedures for structures subjected to fatigue 
loading can be formulated which allows for finite life determination In the 
Initial design stage with reasonable statistical certalnlty. 


DESIGN OF COMPLEX STRUCTURES 

In recent years, finite element stress analysis of complex structures 
subjected to thermal and mechanical loading has reach a high degree of sophis- 
tication and reliability. Now computers perform calculations In mere seconds 
and minutes for problems which only a short time ago would take hours or be 
Impossible to perform. An example of a turbine blade analyzed using finite 
element analysis Is shown In figure 10. Each element of material has asso- 
ciated with It a stress. For design purposes, this stress Is within certain 
established stress or strain limits. The ultimate life of the component Is 
usually based upon empirical calculations or extrapolation from field exper- 
ience. The results are at best highly speculative, by subjecting the com- 
ponent to expensive product Improvement programs (PIP) and by "make and break 
techniques, component lives over a products life time can usually be extended 
to the useful life of the product. 

The key to cost effective design Is to be able to maxlmumlze component 
life and minimize cost at the completion of the final design stage without the 
need for extensive testing or field experience. It Is proposed that this 
objective can be accomplished as follows: 

(1) Having determined the stress of each elemental volume of a component 
using finite-element analysis, It Is possible to establish a life of the 
elemental stressed volume using either equation (3) or (7). 

(2) Using equation (12) It Is then possible to establish the life of the 
entire component. 

(3) Carrying the process a step further, applying the principle of 
equation (8) a "dynamic" capacity of the component can be determined. 



(4) Once the dynamic capacity Is determined, the mission life of a com- 
ponent or system can be established using equation (13). 

(5) Using renewal theory (eqs. (18) and (22)), the number of replacements 
required over a systems useful life and the mean time between failures (MTBF) 
can be determined. 

By design Iteration, design life can be maximized and cost minimized. 
Further, by knowing failure rates, production guidelines can be set prior to a 
product entlrlng the market place. 

In the area of bearing technology, vast technological Improvements 
occurred between Initial publication of Lundberg-Palmgren (1947) and contem- 
porary development. These technological Improvements were factored Into the 
analysis by the use of life adjustment factors (ref. 32) wherein the adjusted 
life L/\ at a 90-percent probabi lity of failure Is 

L a - (D)(E)(F)(G)(H) (£) (23) 

where the life adjustment factors are 

D Materials Factor 

E Processing Factor 

F Lubrication Factor 

G Speed Effect 

H Misalignment Factor 

A similar approach In the design of power-transmitting shafts was taken 
by Loewenthal (ref. 33). The Loewenthal approach Is based on a fatigue limit- 
ing stress based on combined torsion, bending, and axial loading. The fatigue 
limiting stress Is modified by ten factors which In principle are In part 
Incorporated In the Welbull and Lundberg-Palmgren analyses. Those factors 
outside the analysis can be based upon experimental coupon fatigue data. 

Most machine elements are subjected to their maximum design load only for 
a very small faction of their mission life cycle. Hence, by designing to a 
life criterion and a mean cubic load rather than a fatigue limiting stress, It 
Is possible to reduce the size of the machine element without a reduction In 
system reliability. 


CONCLUSION 

A generalized methodology to structural life prediction, design, and 
reliability based upon a fatigue criterion Is advanced. The life prediction 
methodology Is based In part on work of W. Welbull and G. Lundberg and A. 
Palmgren. The approach Incorporates the computed life of elemental stress 
volumes of a complex machine element to predict system life. The results of 
coupon fatigue testing can be Incorporated Into the analysis allowing for life 
prediction and component or structural renewal rates with reasonable statis- 
tical certainty. 
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TABLE 1. - LUNDBERG-PALMGREN FATIGUE LIFE ANALYSIS 


[Based upon 1939 Weibull theory.] 


Analysis 

Year published 

Bearings 

1947 

Spur gears 

1975 

Helical gears 

1976 

Restoraa and refurbished bearings 

1977 

Loaa-life relation 

197b 

Toroidal drive (traction) 

1976 

Optimization of multiroller traction drive 

1980 

Simplified life analysis for traction drives 

1980 

Nasvytis orive 

1981 

Spiral bevel gears 

1982 

Planetary assemblies 

1982 

Transmission assemblies 

1983 

Rotary beam, torsion x tension-compression 

1984 



PERCENT OF BEARINGS FAILED 



Figure 1. - Weibull distribution of bearing fatigue life for 
140-mm bore-size angular-contact ball bearing. 
MIL-L-7808 lubricant; thrust load, 9500 lbs; speed, 

10 000 rpm; temperature, 121 °C (250 °F). 



SHEAR BOND STRENGTH, Ib/sq In. 

Figure 2. - Weibull distribution 
for bond-shear-strength for 
coating of spray powder on 
stainless steel (ref. 11). 
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Figure 3. - Effect oi projected test area on measured and calcu- 
lated bond shear strength of coatings at two surface rough- 
nesses (ref. 11). 
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Figure 4. - Analytically predicted system and component fatigue 
mission lives. Combined material and material processing life 
adjustment factors used for bearings. Gear material is baseline 
material-no factors needed. No lubrication life adjustment 
factors used (ref. 28). 
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Figure 5. - Turboprop reduction gearbox fatigue lives (field 
experience) compared with analytically predicted system 
mission lives. (Material factor implies combined material 
and material processing factor.) (Ref. 28.) 



Figure 6. - Renewal density functions compared for 
rolling-elemont bearings. Also shown is proba- 
bility density for failure of new bearing assuming 
no replacements. Area under curves represents 
probability of failure. 



Figure 7. - Cumulative renewal and failure for rolling- 
element bearings. 
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(b) Soft annealed steel. 


Figure 8. - Comparison between experimental and 
predicted stress-life relation for AISI 52100 steel 
rotating beam specimens for calculated values of 
fatigue limit o u . (Ref. 31. ) 
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(a) Flat beam configuration. 



(b) Stress-life relationship for calculated values 
of fatigue limit o u . 

Figure 9. - Comparison between experimental 
and predicted stress-life relation for flat beam 
specimens made from bearing steel and tested 
in reverse bending (ref. 31). 
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Figure 10. - Finite element analysis 
of turbine engine fan blade. 
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